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ABSTRACT 



Context. Astrophysical jets from various sources seem to be stratified, with a fast inner jet and a slower outer jet. As it is likely that 
the launching mechanism for each component is different, their interface will develop differential rotation, while the outer jet radius 
represents a second interface where disruptions may occur 

Aims. We explore the stability of stratified, rotating, relativistic two-component jets, in turn embedded in static interstellar medium. 
Methods. In a grid-adaptive relativistic hydrodynamic simulation with the AMRVAC code, the non-linear azimuthal stability of two- 
component relativistic jets is investigated. We simulate until multiple inner jet rotations have been completed. 

Results. We find evidence for the development of an extended shear flow layer between the two jet components, resulting from the 
growth of a body mode in the inner jet, Kelvin-Helmholtz surface modes at their original interface, and their nonlinear interaction. 
Both wave modes are excited by acoustic waves which are reflected between the symmetry axis and the interface of the two jet 
components. Their interaction induces the growth of near stationary, counterrotating vortices at the outer edge of the shear flow layer 
The presence of a heavy external jet allows to slow down their further development, and maintain a coUimated flow. At the outer jet 
boundary, small-scale Rayleigh-Taylor instabilities develop, without disrupting the jet configuration. 

Conclusions. We demonstrate that the cross-section of two-component relativistic jets, with a heavy, cold outer jet, is non-linearly 
stable. 

Key words. Stars: winds, outflows - ISM: jets and outflows - Galaxies: jets - methods: numerical, relativity 



1. Introduction 

Several astrophysical jet observations invoke a structured rel- 
ativistic jet, in a direction perpendicular to the jet axis. 
This is the case for Act ive Galactic Nuclei (AGN) jets 
jSiemipinowska et alj I2007h . with sometimes clear evidence 
of a very fast , light jet and a heavy slow external outflow 
dGiroletti et al.l |2004|) . Transversely structured, ultrarelativis- 
tic jet-like outflow is also proposed for Gamma Ray Bursts 
jRossi et al. 2002; Peng et al. 2005), to explain the break in the 
afterglow light curve. 

Theoretical studies for a jet with two componen ts can be 
found in variou s contexts, e.g. for the ca se of AGNinlSoT et alj 
(11989^ and in iRenaud & Heiiril (Il998h . where an electron- 
positron central wind component is surrounded by an exter- 
nal, ideal MHD disk-driven jet. Another two-component o utflow 
model has been proposed by Bogovalov & TsinganosI (l2005h 
where a relativistic pulsar wind is surrounded and collimated 
by an ideal MHD disk-driven wind. Non-relativistic MHD sce- 
narios with two-component jets are also suggested for classical 
T-Tauri stars. In this case of Young Stell ar Objects (YSO), two - 
component models were simulated in ' Meliani eTa D (I2006bl) . 
There, the inner component extracts its energy from a hot 
corona around the central region, while the second component 
is launched from the thin accretion disk. Recent numerical sim- 
ulations have progressed to general rel a tivistic MH D (GRMHD) 
j et lau nching, as done by McKinnevI (12006') and Hard ee et alj 
(12001) . and suggest a similar behavior with two component jets. 
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iMeliani etaP (l2006ah have shown that a GRMHD jet launched 
from the inner region could be collimated by the pressure of ex- 
ternal outflow from the disk. 

The models of disk-driven jets show that while the mass loss 
from the accretion disk is higher than the wind mass loss from 
the inner part, the jet speed increases near the axis, as all the 
models find that the asymptotic velocity of the jet is close to 
the escape speed from the central engine. Indeed, there is a di- 
rect relation between t he asymptotic speed and the depth of th e 
gravitational potential ( Mirabel & Rodrfguezll999l:lLiviol 19991) . 
Another piece of evidence from the models of jet launching that 
the inner component of the jet is faster, is that the disk-driven 
jet mechanism is mainly an ideal MHD mechanism, while the 
wind from the inner region could be more turbulent and is aided 
by non-ideal MHD processes ( Meliani et al.ll200 &5). hence the 
inner jet can be accelerated to high Lorentz facto r, even with 
a weak opening angle of the jet dSpruit et alj|2001h . Therefore, 
a study of the transverse structure of two-component relativis- 
tic jets is a crucial step to understand the physics of high-speed 
jets and the mechanisms responsible in jet launching and colli- 
mation. In fact, the interactions between the two jet components 
may change substantially the dynamics and even ques tion the 
stability of the relativistic jet. Recent stability studies dHarded 
2007) using the linearized RMHD equations, can aid in identi- 
fying the normal modes for magnetized spine-sheath relativistic 
jets. 

While most of these mentioned studies concentrate on mag- 
netized jets, we here intentionally revisit non-linear stability is- 
sues for two-component, relativistic hydro jets. The develop- 
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ment of Kelvin-Helmholtz (KH) instabilities in 3D relativis- 
tic hydro jets has been investigated extensively ('Hardeel l2000t 
[Perucho et al. 2006). Here, we study with unprecedented detail 
the effect of a transverse structure, consisting of an inner, fast, 
and faster rotating outflow, and an external, 'slow' (but relativis- 
tic), heavy, and slowly rotating jet, in turn embedded in a static 
medium. This is done by means of a grid-adaptive computation 
focused on the transverse stability of a cross section of a jet with 
two components. We analyse the effect of density and rotation 
velocity changes at the interface between the two components 
in the jet morphology, and at the outer boundary with the inter- 
stellar medium. As the mechanism responsible for the launching 
of each component is believed to be diff'erent, the energy flux 
and the angular momentum carried in each component is differ- 
ent. We will quantify the nonlinear properties of the structured 
jet and show the ability of the heavy slow jet to collimate and 
stabilize the inner jet. To capture all details, even with a grid- 
adaptive simulation, we perform a 2|D hydrodynamic model, 
as a first step towards more complex models in 3D, eventually to 
be augmented with magnetic fields. Our relativistic hydro simu- 
lation uses a Synge-type equation of state, taking full account of 
effectively varying polytropic index. 



2. Two-component jet model 

We set up a two-component structured jet with a typical total ki- 
netic luminosity flux set to Ljet.Kin = 10"*^ergs/s. The external ra- 
dius of the two-component jet is taken to be /?out ~ O.lpc, a value 
in accord with a jet opening angle 6° at Ipc from the jet source , 
which represent observed values for M87 jBiretta et al.ll2002l) . 
We set the inner jet radius Ri„ - /?out/3, and we ensure that this 
inner component carries only a small fraction /;„ = 0.1% of the 
total kinetic luminosity flux, while the external jet carries the 
remaining /out - 99.9%. For the initial condition, we adopt a 
uniform outflow velocity y, along the jet axis in each jet com- 
ponent, with v'z.out = 0.9428 (hence yz.out ~ 3, since the speeds 
are normalised to the light speed) for the outer, slow jet within 
Rin < R < Rout- This represents a reasonable value for relativis- 
tic jets at the pc scale. The inner jet has a fast outflow set to 
^z.in = 0.9982, with corresponding Lorentz factor -y^.m ~ 16.6. 

Two-component jet models also indicate that the spin of the 
inner beam is higher than the spin of the external jet. Therefore, 
the initial rotation adopts two different profiles, one for the in- 
ner, and one for the external jet. This is inevitable when these 
two jet components are launched from different regions and 
with different mechanisms. We choose for the external jet a ra- 
dially decreasing rotation profile, with a decrease faster than 
Keplerian, as this component is launched from an accretion disk 
and the external streamlines in the jet expand faster than the in- 
ner streamlines. Moreover, we assume that the angular momen- 
tum extracted along each streamline varies little, since we adopt 
a toroidal velocity varying as l/R. For the initial inner jet rota- 
tion profile, we take a toroidal velocity increasing with from 
the axis. This toroidal velocity vanishes on axis and the equiva- 
lent classical centrifugal force given by this profile is constant. 
Summarizing, we take 
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to the small expansion of the inner streamlines, since this com- 
ponent is confined by the external outflow. Also, a fast rotating 
inner jet can effectively extract the angular momentum from the 
central region and carry it away with a very small mass flux. 
This rotation profile of the inner jet obeys the necessary (non- 
relativistic) Rayleigh criterion for stability r{r v^)' > 0, which 
makes the inner jet centrifugally stable. The external outflow is 
marginally stable as r{r v^)' = 0. However, the interface be- 
tween the two components does not verify the Rayleigh crite- 
rion, so the shear flow interface is unstable. 

The densities are also assumed to be constant throughout in- 
ner jet beam, outer jet beam, and external region. The external 
medium is hot and rarified and has a small reference number 
density. The corresponding proper medium density p,ned is in 
the computation used as a scaling value, together with c = 1 
and a unit length of 1 pc. At the parsec scale, the jet head 
has previously shocked the jet beam surroundings, so that the 
jet environment represents a relativistically hot, rarified, static 
external medium. However, in this model we assume a num- 
ber density for the external medium of Icm"^, which is some- 
what higher than the values obtained in some simulations for 
the cocoon surrounding the jet. However, the density of the ex- 
ternal medium has little to no effect in the coUimation of the 
jet cross section as the external jet has slow rotation and us- 
ing smaller value for density in this simulation will only de- 
crease the time step for integration. Values for the jet com- 
ponent proper densities are estimated from the following ar- 
guments prevailing at the jet head. The propagation speed of 
the jet head is estimated using momentum b alance in a ID ap - 



proximation, v. 



,head 
jet 



such that we have a discontinuity in the toroidal velocity at the 
boundary between the two jets when fixing V^ in - 0.05 and 
V^,out = 0.005. A faster rotating inner component could be due 



+ 1) (iMartiet al.l ll997h. 

In this expression, t]r = routPout/iout/Pmed/imed is the ratio be- 
tween the inertia of the (outer) jet and the external medium, 
and hout is the enthalpy of the outer beam. Requiring that 
this outer jet carries a kinetic luminosity flux of /out^'jet.Kin = 
(Tout/iout - I) PoutYoutTT {Rlut - R^u)^'out, we can deduce a relation 
for /lout and thus find the density pout, by prescribing a jet head 
Lorentz factor. We assume that in the interaction with the exter- 
nal medium, the heavy, outer jet undergoes a deceleration to a 
Lorentz factor y^f^ ~ 2.5. This value agrees with the fact that 

' jet.out 

the jet, at the pc scale, is still relativistic, and that the outer jet 
interacts weakly with the external medium. The resulting den- 
sity is Pout ~ 270pmed- A similar argument for the inner, light, 
jet, which therefore must undergo a stronger deceleration, uses 
a Lorentz factor y^^i^^ ~ 1.5, and the density works out to be 

Pin ~ 0.2pn,ed. 

Finally, the pressure is deduced by supposing initial trans- 
verse equilibrium between the pressure gradient and the cen- 
trifugal force, by making the assumption of a constant (but dif- 
ferent for inner versus outer jet) local 'effective polytropic in- 
dex'. Denoting the latter with Tin - 4/3 and Tout = 5/3 (the 
Synge equation of state is not used to determine this initial equi- 
librium), this works out to be 



for R<Rin, 

i \w^, (2) 
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Fig. 1. The effective polytropic index in the initial conditions in Fig. 2. Density in a ID cut through the structured jet after 5 full 
a ID cut through the structured jet. rotations of the inner jet . 
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We choose an initial pressure at the jet axis of po - lO^VmedC 
The governing equ ation of st ate is taken as a Synge-type re- 
lation, also used in Meliani et al.l (.2004). As a result of the above 
profile prescriptions, the external jet is relativistically cold with 
an effective polytropic index Feff = 5/3 and the inner jet has its 
effective polytropic index stratified, since the pressure increases 
fast from the axis where Feff = 5/3, to the inner jet boundary 
where Feff = 4/3. In Fig. [1] we show the effective polytropic 
index of the initial condition. With the ratio between the inertia 
of the outer jet component and the outer shell of the inner jet, 
i.e. youtPout/iout/yfnPin/iin of about 1/2, the local sound speed c, 
is very high in the inner jet and this makes this shell subject to a 
strong interaction with the outer jet. Moreover, with this model, 
the angular momentum carried by the inner jet is of the same 
order as the angular momentum carried by the external jet. The 
outer shell of the inner jet is subsonically rotating with a speed of 
order ~ 0. Ic,, while the outer, faster rotating jet is actually rotat- 
ing at ~ 0.03c5, due to the different thermodynamic conditions. 
Such configuration is favourable to asymmetric instabilities and 
hence their growth times will be of order of the radial sound- 
crossing time (Hardee 2004). 

The computational domain of this simulation is a 2D box of 
size -0.2pc < X < 0.2pc and -0.2pc <y < 0.2pc, and the simu- 
lation is performed in Cartesian c oordinates with the HLLC flux 
formula ("Mi gnone & Bodol2005b . The simulation is ran till time 
t - 22.75. Since the time for a total In rotation of the external 
boundary of the inner jet is 4.4, the simulation follows about 5 
full rotations of the inner jet. The corresponding distance of jet 
propagation of the jet beam during this time is 2 1 .45pc. The sim- 
ulation is done using the AMRVAC code with Synge equation 
of state, using the hybrid block grid adaptivity approach from 
Ivan der Hoist & KeppensI (l2007l) . We take a base resolution of 
128 X 128, allowing for 5 grid levels, reaching effective resolu- 
tion of 20481 



3. Nonlinear evolution and stability 

We first did a ID run, only evolving a radial cut along the cylin- 
drical radius through the jet. This ID simulation demonstrated 



that the two-component jet remains collimated and keeps its in- 
ternal two-component structure. Because of the initial disconti- 
nuity between the two jet components, acoustic waves are cre- 
ated there which propagate through the two jet components. The 
acoustic waves propagating in the inner jet meet at the axis and 
are reflected. The waves propagating in the external jet get par- 
tially reflected on the external boundary of the jet. However, 
these acoustic waves do not disturb the two-component jet struc- 
ture. After three full rotations, the inner jet settles into a station- 
ary state with somewhat higher density on axis, as mediated by 
the propagating acoustic waves. As the outer jet component has 
a high inertia and low sound speed, it efficiently reflects acoustic 
waves at the inner boundary between the two components. 

In 2D, the simulation shows in addition, the development of 
surface and body normal modes. Prior to f ~ 1, the shear flow in- 
terface develops azimuthal Kelvin-Helmholtz instabilities due to 
the velocity gradient. Meanwhile, initially axisymmetric acous- 
tic waves reflect on the jet axis (as in the ID run) and then give 
rise to the growth of a body mode in the inner jet. The latter 
shows a spiral structure with a dominant azimuthal mode num- 
ber m - 4 character Both modes grow and extract energy mainly 
from the inner jet. As mentioned before, the sound speed in the 
external jet is rather slow Cj.om ~ 0.15c, and four times less than 
the sound speed in the inner jet. The body mode in the inner jet 
thus grows fast, as its growth is related to the travel time taken by 
acoustic waves across the inner jet. In such a situation, the body 
mode dominantly disturbs the inner jet component, and the in- 
ner boundary of the jet where the azimuthal surface mode exists. 
This strongly couples the two modes, and from t ~ 1.2 onwards 
the body mode will contribute to surface mode development, and 
induce spiral deformations which penetrate as 4 'arms' into the 
outer jet component. This causes the jet to expand. In this ex- 
pansion phase, the outer surface of the two-component jet be- 
comes unstable to a Rayleigh-Taylor-type instability. This starts 
at about t ~ 3.7 and leads to small-scale structure. However, 
for the entire simulated time-interval, this does not disrupt the 
overall jet structure dramatically. 

The non-linear interaction of the two modes in the inner part 
eventually causes the development of an extended shear flow 
shell, separating the inner and the outer jet. This shear flow 
broadens and compresses the inner jet and thereby increases its 
density, while decelerating it in the vertical (axial) direction due 
to mass conservation. During the formation of the shear flow 
layer, at about t ~ 1 the body mode changes from a dominant 
azimuthal mode number m - 4 character, to an elliptic m - 2. 
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Eventually, as in the ID run, when the acoustic waves in the 
inner jet all but disappear, the inner jet becomes axisymmetric 
again and a more stationary state is reached (at about t ~ 8), 
where the radius of the inner jet drops to Ri^ ~ O.OlSpc while 
the density at the axis reaches pi„ ~ 3.17pmed- From about this 
time, acoustic waves propagating in the extended shear flow re- 
gion get reflected at its bounding interfaces, as the layer becomes 
a kind of resonant cavity. Throughout, a fraction of the angular 
momentum gets transferred to the inner jet region as we find an 
acceleration of the flow rotation to ~ 0.076c at radius Ri„, 
which is needed to radially support the increased pressure of the 
more turbulent shear flow region. This transfer is mostly asso- 
ciated with the 'arms' extending from the rotation axis to the 
outer part of the inner jet, and transfer happens in accord with 
the sound speed. Moreover, the compression of the inner jet 
contributes to the acceleration of the flow rotation. However, the 
simulation shows that there is no transfer of angular momentum 
from the inner jet to the outer jet component, since its rotation 
speed does not increase. All the angular momentum and energy 
extracted from the inner jet ends up in the shear region. 

Fig.[3]shows various quantities at the end time f - 22.75, and 
it can be seen how 4 'arms' have penetrated into the outer jet re- 
gion. The shear flow region then extends between O.OlSpc and 
O.lpc, as seen in Fig. |3] and also in density p in the ID cut 
through the structured jet Fig.|2] The pressure and the centrifu- 
gal force in the sheared shell increases, and together with the ram 
pressure associated with the protruding 'arms', the outer jet re- 
gion has expanded radially. The radius of the outer jet therefore 
rises to Rout ~ O.lSpc. However, the extension of the 'arms' into 
the outer jet region stops when they reach a force equilibrium 
with the external jet, which compresses them. In the process, a 
series of stationary counterrotating vortices develop at the outer 
edge of the shear flow region. These vorticity wave patterns ex- 
tract energy from the inner jet and from waves reflected at both 
bounding interfaces of the formed shear flow region. In the sim- 
ulation, we find that the spiral arms and the vorticity waves are 
coupled, and the vortices remain near stationary. 

The thermal energy in the inner jet remains relativistic, with 
an effective polytropic index Fgff ~ 1.38, while in the outer 
jet thermodynamic conditions remain classical with an effective 
polytropic index Feff ~ 1.652. However, the thermal energy in 
the shear flow layer, and in the vortex structures is mildly rel- 
ativistic. In the shear flow the effective polytropic index is of 
the order Feff ~ 1.45, and in the center of the vortex we find 

Feff ~ 1.5. 
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4. Conclusions 

We investigated the dynamical consequences of an inner fast, 
light, and fast rotating outflow, on a slower, heavy and slower 
rotating relativistic jet. The latter could be launched from an ac- 
cretion disk surrounding the engine for the inner jet. We per- 
formed a numerical simulation to study the azimuthal stability 
of the cross section of this two component jet. The massive, cold 
outer jet makes the inner, light jet act as a resonant cavity due to 
strong reflection of acoustic waves at the interface. This induces 
various linear instabilities which interact strongly. The simula- 
tion shows evidence for a transition to a near-stationary state, 
three-component state, with a separating shear flow region be- 
tween the two components. In the process, stationary counter- 
rotating vortices develop. Their growth gives rise to significant 
thermal and ram pressure in the shear region, which compresses 
the inner jet and partially expands the external jet. However, the 
heavy, cold outer jet slows down the instability development. 
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Fig. 3. Density, toroidal velocity, effective polytropic index, and 
Lorentz factor for the two-component jet, after 5 full rotations of 
the inner jet. 
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and the overall increased cross section remains coUimated after 
several dynamical rotation times of the inner component. The 
resulting more turbulent shear layer between the two jet compo- 
nents could be responsible for continuous in-situ reacceleration 
of electrons, and have a mea surable impact on the radiatio n spec- 
tra, as already suggested bv lStawarz & Ostrowskil (|2002|) . 

The non linear stability analysis of the jet cross section pre- 
sented here allows to investigate precisely the development of 
azimuthal instabilities and their influence on the inner and the 
external structure of the jet cross section. The azimuthal insta- 
bilities between the two jet components ultimately modify the 
speed of the inner jet and the width of the external jet which 
should influence the 3D structure of the jet. In future work, we 
will study the interaction between two components in full 3D 
simulations. They should explore the influence of the initial ro- 
tation profile and the relative influence of azimuthal versus lon- 
gitudinal instabilities for realistic multi-component jets. Finally, 
magnetic fields must be included for quantifying their further 
role in relativistic structured jet collimation. 
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